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APPLICATION  OF  THE  PUGH.  EICHELBERG  AND  ROSTOKER 


THEORY  TO  THE  MRL  38  MM  SHAPED  CHARGE 


1 .  INTRODUCTION 


The  use  and  understanding  of  shaped  charge  waheads  has  been,  and 
continues  to  be,  of  considerable  Interest  to  MRL.  The  Australian  Amy 
currently  employs  two  shaped  charge  weapons,  the  66  mm  M72  L1A2  (LAN)  and  84 
am  Carl  Gustaf  anti-tank  weapons,  and  there  Is  continuing  Interest  In  the  103 
mm  MILAN,  a  shaped  charge  warhead  irtilch  Incorporates  wave  shaping  for 
Increased  performance.  The  selection,  performance,  and  possible  modification 
of  such  weapons  to  fit  Australian  requirements  is  obviously  facilitated  by  a 
thorough  understanding  of  the  mechanisms  of  liner  collapse  and  Jet  formation. 


A  38  mm  shaped  charge  with  42*  angle  copper  liner  and  Comp.  B 
explosive  was  designed  and  tested  at  MRL  some  years  ago  [1].  This  particular 
shaped  charge,  now  known  as  the  MRL  standard  shaped  charge.  Is  currently  being 
used  by  M.C.  Chick  and  co-workers  for  a  very  Illuminating  Investigation  on  the 
Jet  Initiation  and  disruption  of  bare  and  covered  Comp.  B  and  other  explosives 
[2a, 2b]. 


The  process  of  liner  collapse  and  Jet  formation  In  the  MRL  Standard 
Shaped  Charge  has  recently  been  modelled  by  D.L.  Smith  using  the  two- 
dimensional  hydrodynamics  finite-difference  code  HELP  [3] ,  and  Initial 
comparison  between  flash  radiographs  of  liner  collapse  and  the  HELP  output  are 
most  encouraging  [4],  although  the  code  does  have  problems  In  accurately 
modelling  the  Jet  tip.  While  HELP  appears  to  be  capable  of  modelling  the  Jet 
formation  process  for  the  MRL  standard  shaped  charge  with  good  overall 
precision,  it  Is  not  Ideally  suited  to  parametric  or  design  studies  because  of 
the  extensive  amount  of  computing  time  required,  with  each  run  costing  several 
hundred  dollars. 


A  different  approach  to  the  analysis  of  shaped  charge  performance  Is 
to  use  one  of  the  many  so  called  **one  dimensional”  codes  (ID),  all  of  which 
are  based  on  the  analytical  work  of  Pugh,  Elchelberg  and  Rostoker  [ 5] , 

(PER).  A  number  of  these  codes  have  appeared  In  the  last  few  years,  the  best 


known  of  these  probably  being  the  BASC  code  of  J.T.  Harrison  of  BRL  [6],  while 
others  are  the  OESC-1  code  of  Car leone  et  al  of  Oyna  East  Corp*  [7],  the 
JETFORM  code  at  RARDE  [8],  and  the  TB/ISL  code  of  Hennequin  in  France  [9]. 

Each  of  these  codes  is  based  on  the  evaluation  of  siiiq>le  analytical 
expressions  rather  than  the  solution  of  a  set  of  one  dimensional  finite 
difference  equations,  as  the  name  would  suggest,  and  hence  the  saving  in 
computer  time  when  compared  with  a  2D  hydrocode  run  is  enormous.  Such  codes 
are  ideally  suited  to  parametric  studies  and  are  an  integral  part  of  the 
overall  design  process,  as  recently  outlined  by  Carleone  et  al  [10]. 

Despite  having  acquired  several  large  2D  hydrocodes  in  the  past  few 
years,  the  most  notable  of  these  being  the  previously  mentioned  HELP  code,  as 
well  as  2DL  [11],  HULL  [12],  and  HEMP  [13],  MRL  has  not  yet  obtained  any  of 
the  ID  codes  for  shaped  charge  analysis.  It  was  felt  that  an  MRL  ID  code 
should  be  developed  first  and  applied  to  the  standard  MRL  shaped  charge  before 
any  of  the  overseas  ID  codes  were  considered.  The  advantages  of  such  an 
approach  are  an  Increased  understanding  of  the  basic  physics  Involved  in  the 
liner  collapse  and  jet  formation  process,  and  a  greater  appreciation  of  what 
to  look  for  in  other  ID  codes  if  these  become  available  at  MRL.  The  purpose 
of  this  report  then  is  to  describe  the  progress  made  so  far  in  the  modelling 
of  the  MRL  standard  shaped  charge  using  the  **1D  code**  type  of  approach.  In 
the  next  section  a  description  is  given  of  the  basic  equations.  No  attempt 
has  been  made  to  derive  these  as  a  detailed  description  of  the  derivation  of 
the  basic  PER  equations  has  recently  been  given  in  a  companion  MRL  report 
[14],  which  provides  the  essential  background  for  the  understanding  of  the 
work  described  here.  The  modelling  of  a  typical  shaped  charge  using  the 
approach  described  in  this  report  requires  the  determination  of  three 
phenomenological  constants,  which  are  found  by  fitting  to  some  experimental 
data  for  the  BRL  105  mm  unconflned  shaped  charge.  This  process  is  described 
in  Section  3.  In  Section  4  we  discuss  any  changes  which  have  to  be  made  to 
these  constants  to  enable  us  to  model  related  shaped  charge  systems,  and  we 
then  apply  the  model  to  the  MRL  38  mm  shaped  charge.  Predictions  of  the  model 
are  compared  with  the  recent  experimental  results  of  M.C.  Chick  [15]  and  the 
flash  radiographs  of  D.L.  Smith  and  l.B.  MacIntyre  [4].  The  effect  of  varying 
both  liner  thickness  and  cone  angle  on  the  velocity  of  the  jet  tip  is  then 
Investigated.  Finally,  in  Section  5,  we  discuss  additional  possible 
applications  of  the  work  described  here,  and  several  methods  for  refining  the 
model. 


2.  DESCRIPTION  OF  EQUATIONS 

Before  the  PER  equations  can  be  used  to  predict  the  velocities  of 
jet  and  slug  an  expression  must  be  found  for  the  velocity  Imparted  to  the 
liner  by  the  explosive  loading.  This  velocity  is  characterized  by  its  speed 
Vq  and  the  angle  fi  between  the  velocity  vector  and  the  perpendicular  to  the 
original  liner  surface  at  each  point  along  its  length.  These  two  variables 
are  related  by  the  Taylor  equation  [20],  which  for  steady  state  conditions  is 


where  U  >*  Uq/cos  a,  Uq  Is  the  velocity  of  detonation,  and  a  Is  the  cone  half 
angle.  Equation  (1)  assumes  steady  state  detonation  In  the  direction  of  the 
charge  axis.  It  has  been  generalized  to  non  steady  conditions  by  Chou  et  al 
[16],  but  we  shall  not  be  concerned  with  this  here. 

As  described  In  a  previous  report  [14],  there  are  basically  two 
approaches  to  the  determination  of  We  can  calculate  the  bending  angle 

5  using  a  form  of  the  Richter  equation  [17]  and  then  use  the  Taylor  equation 
to  find  V^,  or  we  can  first  calculate  using  an  advanced  Gurney  equation 
[18]  and  then  find  5  from  equation  (1).  Both  methods  were  tried  for  the  MRL 
38  mm  shaped  charge  and  the  first  method  was  found  to  be  the  simplest  In 
practice  for  the  accuracy  required.  6  is  calculated  using  an  equation  of  the 
form 


^  +  Kp  (  £  )  ,  (2) 

^O 


where  p  and  e  are  the  density  and  thickness  of  the  liner  material  and  e  Is  the 
thickness  of  the  explosive.  ^  and  K  are  phenomenological  constants  which 
are  found  by  fitting  to  a  well  characterized  BRL  105  mm  shaped  charge.  This 
is  described  In  the  next  section.  Equation  (2)  is  the  Richter  equation  as 
modified  by  Defourneaux,  and  more  Information  on  this  can  be  found  In 
reference  [19].  Harrison  [6]  also  has  modified  equation  (2)  to  take  Into 
account  the  effect  of  different  case  confinements.  As  the  BRL  105  mm  is 
unconflned  and  the  MRL  38  mm  only  lightly  confined  It  Is  assumed  here  that 
this  correction  can  be  neglected. 

Once  V  and  6  have  been  calculated  the  PER  equations  can  be  used  to 
calculate  the  collapse  angle  8  and  the  velocities  and  masses  of  Jet  and 
slug.  To  calculate  8  a  modification  of  the  original  PER  equation  Is  used. 

The  modified  expression  Is 

tan  8  *  tan(8'*^  +  A8),  (3) 


where  68  is  defined  by 


tan  A8  “  -  (xslna/cos<a+5 )cos6)  (V^j'/V^)  (4) 

8^  is  the  collapse  angle  which  Is  applicable  In  the  steady  state  case  and  was 
first  calculated  by  Blrkhoff,  MacDougall,  Pugh  and  Taylor  [20].  This  Is 
given  by 


X  is  a  measure  of  distance  along  the  cone  axis  from  the  apex  and  the  prime 
denotes  differentiation  with  respect  to  x.  Equation  (3)  was  recently  derived 
by  Hlrsch  [21]. 

Once  0  has  been  calculated  the  velocities  of  each  element  of  Jet  and 
slug  are  calculated  from 

V.  -  V  coa(a  +  «  -  0/2)  /  sin(0/2)  (6) 

j  o 

V  -  V  aln(a  +  6  -  0/2)  /  cos(0/2)  ,  (7) 

s  o 

and  the  masses  from 

mj/m  ”  8in^(0/2)  (8) 

m^/m  -  co8^(0/2),  (9) 

where  m  is  the  liner  mass  per  unit  length.  These  equations  were  originally 
derived  by  Pugh,  Elchelberg  and  Rostoker  [5]  and  are  discussed  In  detail  In 
reference  [14] . 

When  modelling  both  the  BRL  105  mm  and  MRL  38  mm  shaped  charges  It 
was  found  necessary  to  take  account  of  the  inverse  velocity  gradient  effect, 
(IVG).  This  was  not  Included  In  the  original  PER  equations  but  has  been 
discussed  In  detail  by  Carleone  et  al  [22].  Here  the  effect  Is  modelled  by 
assuming  that  the  velocity  of  each  liner  element  after  being  hit  by  the 
detonation  wave  has  the  following  form 

V  -  (1  -  exp(-(t  -  T)/t),  (10) 

where  V  Is  the  velocity  at  any  Instant,  Is  the  final  velocity,  T  Is  the 
time  at  which  the  detonation  front  reaches  the  element,  and  t  Is  a  time 
constant  (yet  to  be  determined)  for  the  acceleration.  Whilst  equation  (10) 
is  obviously  Incorrect  for  very  small  times.  Its  overall  applicability  for 
liner  motion  has  been  amply  justified  by  Chou  et  al  [16]. 

The  calculations  described  In  this  report  assume  that  the  velocity 
of  the  jet  tip  has  already  been  experimentally  determined.  This  quantity, 
together  with  the  velocity  of  detonation  Uq  and  the  cone  half  angle  a,  are 
sufficient  to  give  Initial  estimates  of  V^,  6  and  0  using  the  original  steady 
state  equations  of  Blrkhoff  et  al,  l.e.  equations  (1),  (5)  and  (6).  These 
equations  can  be  solved  to  give  the  following  expression  for  V 


tan  (a/2) 

(  tan  (o/2)  +  1)  ) 


O' 


Equation  (1^  Is  then  used  to  calculate  d,  and  then  equation  (5)  used  to 
calculate  g  . 

The  calculation  starts  by  dividing  the  shaped  charge  liner  Into  NN 
segnents,  where  NN  Is  typically  30  or  40.  Hie  bending  angle  for  each  element 
Is  calculated  from  equation  (2),  Is  calculated  from  equation  (1),  and  then 
the  collapse  angle  E  Is  calculated  from  equation  (3).  At  this  stage  various 
other  quantities  of  Interest  are  also  calculated.  These  Include  T*,  which  Is 
the  time  It  takes  an  element  of  the  liner  to  reach  the  axis  after  being  hit  by 
the  detonation  wave,  z,  which  Is  the  point  along  the  liner  axis  at  which  each 
element  collides  to  form  Jet  and  slug,  and  DIST,  which  Is  the  distance  each 
element  travels  before  hitting  the  axis.  z  and  DIST  are  given  by 


-r. 

-■s-V-'n  .vv 
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z  ■  x(l  +  tana  tan(a  +6))  , 


DIST  -  x  taaa/cos(a  +6)  . 


These  expressions  were  first  given  by  Carleone  and  Chou  [23] . 

The  IVG  effect  Is  now  Included  by  assuming  that  each  eleir^fit  has  a 
velocity  given  by  equation  (10).  The  position  as  a  function  of  time  for  each 
element  Is  now  calculated  until  the  distance  travelled  equals  DIST.  At  this 
stage  T*  Is  reset  to  the  current  time  and  V  Is  reset  to  the  current  V.  A 

new  6  Is  calculated  from  equation  (1)  and  then  a  new  S  from  equation  (3). 

Is  then  calculated  from  equation  (6)  and  slower  leading  Jet  elements  are 
allowed  to  collide  assuming  perfectly  plastic  Impacts.  The  equilibrated  Jet 
tip  velocity  is  then  calculated  from  conservation  of  mass  and  momentum. 


In  the  next  section  we  describe  the  application  of  this  approach  to 
the  BRL  105  mm  unconflned  shaped  charge. 


3.  DETERMINATION  OF  PHENOMENOLOGICAL  CONSTANTS 


Before  using  the  scheme  outlined  In  the  previous  section  to  model 
the  MRL  38  mm  shaped  charge  we  have  to  determine  the  constants,  R,  4  and  t. 
As  In  any  phenomenological  theory  these  are  found  by  fitting  to  experimental 
data.  Because  the  BRL  105  mm  unconflned  shaped  charge  has  been  studied  so 
extensively,  and  because  of  Its  apparent  scaling  to  the  MRL  shaped  charge.  It 
was  decided  to  use  the  available  data  on  this  charge  to  determine  the 
constants. 


The  dimensions  of  the  charge  are  depicted  in  Fig.  1.  For  modelling 
purposes  the  true  charge  geometry  is  replaced  by  the  pointed  cone  shown  by  the 
dashed  line  in  Fig.  1.  The  radius  is  taken  to  be  43.2  mm  and  the  thickness  of 
the  copper  liner  to  be  2.7  mm.  The  velocity  of  detonation  in  the  Comp.  B  is 
7.98  mm/ps  and  the  experimental  value  for  the  velocity  of  the  jet  tip  is 
7.0  mm/ps.  All  data  have  been  taken  from  [6].  Fig.  2  shows  the  experimental 
values  for  and  V.  as  a  function  of  scaled  length  of  the  cone.  Also  shown 
is  the  theoretical  tit  to  this  data  using  the  BRL  one-dimensional  BASC  code  of 
Harrison  [6] .  The  experimental  points  are  from  the  work  of  Allison  and 
Vltall  [24].  Fig.  3  shows  the  experimental  values  for  cumulative  mass  of  Jet 
versus  cumulative  mass  of  cone  and  also  the  BASC  code  prediction.  As  can  be 
seen,  even  though  the  shapes  are  similar,  there  is  considerable  discrepancy 
between  theory  and  experiment. 

To  model  this  shaped  charge  using  the  ideas  outlined  in  the  previous 
section  we  begin  by  arbitrarily  choosing  values  for  K  and  ^  and  then 
comparing  the  resulting  V,  and  values  with  the  experimenBal  results.  After 
a  period  of  trial  and  error  the  values  K  •  1.5  and  ^  «  28*  were  decided 

upon.  These  give  a  reasonable  fit  to  the  V.  and  Sata,  as  shown  in  Fig.  4, 
and  also  quite  good  agreement  with  the  BASC^result  for  the  cumulative  mass  of 
jet  versus  cumulative  mass  of  cone,  as  shown  in  Fig.  3  by  the  solid  curve 
labelled  t  >■  0. 

The  next  step  is  to  include  the  IVG  effect.  This  is  done  by 
choosing  T  to  have  the  value  1.0  and  then  slowly  increasing  this  value  until 
better  agreement  is  obtained  with  the  experimental  results.  Fig.  5  shows  the 
effect  on  when  t  *  2.5.  The  IVG  effect  has  reduced  the  value  of  V  near 
the  apex  of  the  cone  so  that  it  now  agrees  with  the  BASC  result,  although 
there  is  still  quite  a  discrepancy  with  experiment.  By  increasing  t  to  8.0 
this  discrepancy  is  reduced  to  an  acceptable  level,  as  can  also  be  seen  from 
Fig.  5.  For  T  “  8.0  the  equilibrated  jet  tip  velocity  is  7.0  mm/ps,  in  good 
agreement  with  the  experimental  results,  and  the  plot  of  cumulative  mass  of 
jet  versus  cumulative  mass  of  cone  also  shows  much  better  agreement  with 
experiment.  This  is  shown  in  Fig.  3. 

The  level  of  agreement  between  theory  and  experiment  shown  here 
suggests  that  the  BRL  105  mm  unconflned  shaped  charge  is  satisfactorily 
modelled  by  the  scheme  outlined  in  the  previous  section  with  the  values  of  the 
constants  just  given.  We  note  also  that  this  level  of  agreement  is 
independent  of  the  number  of  segments  into  which  the  shaped  charge  is 
initially  divided.  In  the  next  section  we  apply  this  approach  to  the  MRL 
38  mm  shaped  charge. 


4.  APPLICATION  TO  MRL  38  MM  SHAPED  CHARGE 


The  dimensions  and  composition  of  the  MRL  shaped  charge  are  shown  in 
Fig.  6.  For  modelling  purposes  we  take  the  radius  of  the  liner  to  be  19.1  mm 
with  a  thickness  of  1.02  mm.  The  velocity  of  detonation  in  the  Australian 


Comp.  B  Is  7.76  mm/ps  [31]  and  the  best  experimental  value  for  the  velocity  of 
the  jet  tip  Is  7.33  mm/ps  [4]. 

Before  proceeding  further  we  now  have  to  examine  how  the  three 
phenomenological  constants  scale  with  the  size  and  type  of  shaped  charge  under 
consideration.  Equation  (2)  was  originally  derived  for  planar  geometries. 
Under  these  conditions  K  and  ^  depend  only  on  the  particular  exploslve/metal 
combination  and  are  Independent  of  size.  Both  Defourneaux  [19]  and  Harrison 
[6]  point  out  that  ^  does  depend  on  the  angle  of  Inclination  of  the 
detonation  wave  to  tRe  metal  surface,  but  as  the  liner  half-angle  Is  21*  for 
both  the  BRL  105  mm  and  MRL  38  mm  shaped  charges  this  Is  of  no  concern  to 
us.  The  application  of  equation  (2)  to  the  curved  geometry  of  shaped  charges 
Is  rather  outside  the  scope  of  Its  Initial  derivation,  but  this  approximation 
has  been  used  now  In  many  papers  on  shaped  charge  modelling  and  the  results 
have  Indicated  Its  applicability  In  this  area.  We  conclude  then  that  the 
values  of  the  constants  K  and  ^  derived  for  the  BRL  105  mm  shaped  charge 
should  also  be  applicable  to  the  MRL  38  mm  shaped  charge,  although  the  smaller 
run-distance  to  the  apex  of  the  cone  may  affect  the  Impulse  delivered  to  the 
liner  and  therefore  have  a  sleight  affect  on  the  constants.  This  Is  currently 
under  Investigation  [27].  The  smaller  run-distance  will  also  mean  that  the 
detonation  front  Is  not  perfectly  planar  and  so  the  angle  of  Inclination  will 
not  be  exactly  21*,  but  the  effect  of  this  on  the  constants  K  and  ^  Is 
expected  to  be  very  small. 

The  question  with  regard  to  t  Is  less  clear.  t  does  In  fact  vary 
with  position  along  the  liner  and  Chou  et  al  [16]  give  expressions  for  this  In 
their  paper.  Our  t  Is  position  Independent  and  Is  an  Integrated  form  of  the 
one  used  by  them,  and  consequently  Its  dependence  on  the  overall  size  of  the 
shaped  charge  Is  harder  to  predict.  For  the  moment  we  simply  regard  t  as  a 
free  parameter  to  be  fixed  by  requiring  that  the  equilibrated  jet  tip  velocity 
be  equal  to  the  experimental  jet  tip  velocity.  This  results  In  a  value 
for  T  of  2.3  us. 

The  predictions  of  this  modelling  scheme  when  applied  to  the  MRL 
shaped  charge  are  now  described.  Fig.  7  shows  the  velocity  gradient  within 
the  jet  before  equilibration  of  the  jet  particles  has  occurred.  The  liner  was 
Initially  divided  Into  30  segments  for  this  calculation,  but  the  result  Is 
Independent  of  the  number  of  segments  chosen.  The  Inverse  velocity  gradient 
effect  Is  clearly  evident,  and  calculations  allowing  fully  plastic  Impact  of 
the  jet  particles  show  that  material  from  the  first  ten  segments  of  the  cone 
collide  to  from  a  compact  jet  tip  particle  with  a  mass  of  0.27  g.  The  data 

contained  In  Fig.  7  are  replotted  In  Fig.  8  after  equilibration  of  the  jet  tip 
velocity  to  show  the  velocity  gradient  within  the  jet  as  a  function  of  the 
scaled  length  of  the  jet.  Also  shown  In  Fig.  8  are  experimental  values  for 
the  velocity  gradient  obtained  from  recent  experiments  by  M.C.  Chick  at  MRL 
[15].  These  were  obtained  by  taking  flash  radiographs  of  the  jet  after  It 
had  broken  up  into  a  number  of  small  segments.  The  velocity  of  each  segment 
was  then  obtained  by  measuring  the  distance  It  had  travelled  between  two  X-ray 
flashes.  Unfortunately  the  X-ray  cassettes  were  not  large  enough  to  measure 
all  of  the  jet  particles  produced  and  so  an  estimate  of  the  total  number  had 
to  be  made.  Several  BRL  reports  (such  as  [25,26])  have  noted  that  the  number 
of  particles  Into  which  a  jet  will  break  Is  Independent  of  the  size  of  the 
shaped  charge  and  depends  only  on  the  explosive,  provided  all  dimensions  are 


scaled  appropriately*  For  a  copper  lined  shaped  charge  with  a  21**  half  angle 
and  Comp*  B  explosive  the  jet  breaks  Into  about  50  particles  having  velocities 
within  the  range  ^ A(t±n)  (approximately  7*5  mm/ps)  down  to  2*0  mm/ps* 

Particles  with  velocities  less  than  2*0  mm/ps  are  also  produced  close  to  the 
back  of  the  jet,  but  these  are  not  counted  by  BRL  as  they  are  Incapable  of 
penetrating  a  target  [25]*  With  the  above  In  mind,  a  plot  of  particle 
velocity  versus  particle  number  was  made  for  the  data  obtained  by  M*C. 

Chick*  This  Is  shown  In  Fig*  9*  A  maximum  number  of  38  jet  particles  could 
be  counted  for  this  particular  set  of  experiments  and  the  velocity  of  the 
slowest  particle  measured  was  approximately  4  mm/ps*  A  straight  line  drawn 
through  the  data  and  then  extrapolated  back  to  a  particle  velocity  of 
2  mm/ps  Indicates  a  particle  number  of  55,  which  Is  In  good  agreement  with  the 
BRL  data  as  the  38  mm  shaped  charge  may  be  too  small  for  exact  scaling  to 
apply*  The  velocity  of  the  tall  end  of  the  jet  predicted  from  Fig*  8  is 
0*18  mm/ps  and  from  Fig*  9  this  corresponds  to  72  particles*  Assuming  this  Is 
the  total  number  of  particles  produced  in  the  breakup  of  the  jet  we  can  then 
plot  the  experimental  fragment  velocities  as  a  function  of  the  scaled  length 
of  the  jet  as  shown  In  Fig*  8*  The  level  of  agreement  between  experiment  and 
the  prediction  of  the  model  Is  quite  pleasing  considering  the  uncertainty 
involved  In  analysing  the  experimental  data  and  the  simplicity  of  the  model* 
Obviously  this  agreement  depends  fairly  strongly  on  the  estimate  of  the  total 
number  of  particles  produced  on  jet  breakup,  and  further  experiments  are 
planned  to  enable  this  number  to  be  determined  much  more  accurately  [27]. 

We  now  consider  the  collapse  angle  8  as  a  function  of  time  (t  >  0 
corresponding  to  the  firing  of  the  detonator)*  Fig*  10  shows  the  prediction 
of  the  model  and  a  comparison  with  the  experimental  results  obtained  by 
O.L*  Smith  and  I*B*  MacIntyre  [4]*  It  should  be  noted  that  the  estimation  of 
these  collapse  angles  from  the  flash  radiographs  Is  a  particularly  difficult 
measurement  to  make,  and  Is  subject  to  appreciable  error  because  of  the 
difficulty  in  estimating  the  position  of  the  inner  edge  of  the  liner*  The 
results  shown  in  Fig*  10  were  obtained  by  D*L*  Smith  and  I*  MacIntyre  [4] 
using  Image  enhancement  techniques  and  are  probably  the  best  which  can  be 
obtained,  but  even  so  the  measurements  of  the  larger  angles  are  subject  to 
such  large  uncertainties  that  they  have  not  been  Included  here*  The 
agreement  between  theory  and  experiment  Is  again  quite  acceptable  considering 
the  simplicity  of  the  model  and  the  experimental  errors* 

The  above  results  indicate  that  the  equations  described  in  Section  2 
provide  a  realistic  model  for  linear  collapse  and  jet  formation  In  the  MRL  38 
mm  shaped  charge,  and  so  we  can  now  use  this  model  to  Investigate  various 
parametric  dependencies*  For  example.  Fig*  11  shows  the  effect  of  varying  the 
cone  half  angle  a  on  jet  tip  velocity*  The  dashed  line  has  been  calculated 
using  the  simple  steady  state  theory  of  Blrkhoff  et  al  [20],  while  the  full 
line  has  been  calculated  using  the  present  model*  The  dots  represent 
experimental  values  obtained  by  BRL  for  their  1-1/2  Inch  diameter  copper  lined 
shaped  charge  with  Comp  B  explosive  [26]*  The  experimental  points  show  that 
as  the  cone  half  angle  a  Increases  then  \  decreases*  This  trend  Is 

predicted  both  by  the  steady  state  theory^and^^the  model  developed  here,  with 
the  latter  giving  better  agreement  with  the  experimental  values  as  a  Increases 
from  10“  to  30“ *  Beyond  this  value  of  a  the  model  predicts  that 
becomes  almost  constant,  while  the  simple  theory  predicts  a  continuous*^ 
decrease  in  y  The  reason  for  this  Is  that  as  a  increases  each  element 

on  the  liner  near'^the  apex  of  the  cone  has  a  greater  distance  to  travel  and  Is 


therefore  able  to  obtain  a  greater  fraction  of  Its  terminal  velocity  before 
colliding  on  the  axis.  This  Increase  In  element  velocity  on  hitting  the  axis 
balances  the  decrease  In  V ..  .  .  due  to  the  changing  geometry  and  so  . 

remains  approximately  constant.  The  experimental  point  at  a  ~  45*  shows  that 
V.^ti  \  does  Indeed  not  decrease  as  much  as  predicted  by  the  simple  theory, 
but  Ir  does  decrease  and  not  remain  constant.  The  reason  for  the  failure  of 
the  model  In  this  region  Is  because  of  the  neglect  of  the  dependence  of  the 
constants  K  and  ^  In  equation  (2)  on  the  angle  of  Inclination  of  the 
detonation  wave  to  the  liner,  1.  Fig.  2  of  reference  [6]  shows  that  as  1 
decreases  (l.e.  a  Increases)  both  K  and  4  become  rapidly  varying  functions  of 
1.  The  model  can  be  fitted  to  the  experimental  point  at  a  ~  45*  simply  by 
changing  ^  from  28*  to  20*,  which  Is  In  agreement  tilth  the  trend  shown  In 
Fig.  2  of  reference  [6]. 

The  model  could  also  be  used  to  Investigate  other  parametric 
dependencies,  such  as  the  variation  of  .  with  the  thickness  of  the 

liner.  Here  It  would  first  be  Important '’to  iind  the  dependence  of  t  on  the 
thickness  of  the  liner,  and  this  could  be  done  either  by  analysing  the  trends 
from  some  of  the  BRL  data,  or  using  a  code  such  as  2DL  [11]  to  explicitly 
calculate  t  for  various  liner  thicknesses  (and  also  for  different 
explosives).  An  application  of  the  model  to  wave  shaping  In  the  MRL  shaped 
charge  has  just  been  completed  and  will  be  reported  elsewhere  [32]. 


5.  CONCLUDING  REMARKS 


A  simple  model  of  liner  collapse  and  jet  formation  for  the  MRL  38  mm 
shaped  charge  has  been  described  and  the  predictions  of  the  model  have  been 
found  to  compare  favourably  with  recent  experimental  results.  The  model  could 
obviously  be  Improved  In  several  ways  If  more  accuracy  was  required,  and  two 
of  these  options  are  now  described: 

(1)  The  calculation  of  the  Impulse  exerted  on  the  liner  by  the 

detonation  wave  could  be  made  using  an  Improved  Gurney  model  of  the  type 
described  by  Chanteret  [18]  or  Hennequln  [28]. 

This  would  have  the  advantage  of  using  known  physical  constants  for 
the  particular  exploslve/llner  combination  and  would  not  require  fitting  to 
phenomenological  constants  such  as  K  and  4  .  An  even  more  sophisticated 
method  of  calculating  liner  Impulse  would  §e  to  use  a  2D  code  such  as  HEMP  or 
2DL  to  provide  the  values  of  V  and  5  to  feed  Into  the  present  model.  Such  an 
approach  has  been  used  by  Van  Thiel  and  Levitan  [29]. 

(11)  The  modelling  of  the  jet  tip  could  be  Improved  and  a  plotting 

program  could  be  written  to  draw  the  outline  of  the  jet,  slug,  and  liner  at 
selected  times.  The  simplest  way  to  do  this  would  be  to  assume  that  the 
accumulated  segments  In  the  jet  tip  form  a  spherical  blob,  but  already  we  know 
that  this  Is  not  a  very  accurate  description  of  the  process.  The  jet  tip 
appears  flattened  along  Its  leading  edge,  and  there  Is  evidence  that  ablation 
has  a  significant  effect  on  the  shape  of  the  jet  tip.  It  would  be  of  Interest 
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to  understanding  this  formation  process  In  more  detail,  however,  because  of 
Its  obvious  relevance  to  the  jet  Initiation  work  of  M.C.  Chick  [2] . 


Other  possible  Improvements  are  the  Incorporation  of  the  non  steady 
Taylor  angle  formulation  of  Chou  et  al  [16],  and  the  Inclusion  of 
compressibility  along  the  lines  described  by  Chou,  Carleone  and  Karpp  [30] . 

In  Its  present  form  the  model  describes  experimental  observations  well  and  has 
demonstrated  good  capabilities  for  parametric  studies,  thus  further  refinement 
would  have  to  be  justified  by  more  detailed  experimental  characterisation  of 
shaped  charges. 
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FIGURE 


Experimental  values  for  Vj  and  for  the  BRL  105  mm  shaped 
charge.  The  solid  llnes''are  the  predictions  of  the  BASC  code 
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CUMULATIVE  MASS  OF  CONE  0 


FIGURE  3  Cumulative  mass  of  Jet  versos  cumulative  mass  of  cone  for  the  BRL 
105  mm  shaped  charge.  The  squares  are  the  experimental  values 
while  the  circles  are  the  predictions  of  the  BASC  code.  The 
solid  line  Is  the  result  predicted  by  the  MRL  model  without 
Inclusion  of  the  inverse  velocity  gradient  effect,  while  the 
dashed  line  Is  the  prediction  when  this  effect  is  Included,  with 
T  ■  8.0. 


•  • 


((ir/uiui  ' 


IP 


SCALED  LENGTH  OF  JET 


Velocity  gradient  within  the  MRL  jet  as  a  function  of  scaled 
length  of  Jet.  The  solid  line  Is  the  prediction  of  the  model, 
the  experimental  points  are  from  recent  work  by  M.C.  Chick  [15] 


PARTICLE  NUMBER 


Particle  velocity  versus  particle  number  for  the  MRL  jet  after 
breakup.  The  different  data  points  represent  two  experiments 
from  recent  work  by  M.C.  Chick  [15]. 


FIGURE  10  Collapse  angle  0  for  the  NRL  shaped  charge  as  a  function  of  time 

after  firing  the  detonator.  The  solid  curve  is  the  prediction  of 
the  model,  the  experimental  points  are  from  recent  work  by 
D.L.  Smith  and  I.B.  MacIntyre  [4].  The  dashed  line  represents 
the  time  taken  by  the  detonation  wave  to  travel  from  the  detonator 
to  the  cone  apex. 
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FIGURE  11  Velocity  of  jet  tip  as  a  function  of  cone  half  angle.  The  dashed 
line  has  been  calculated  from  the  steady  state  theory  of  Blrkhoff 
et  al  [20] ,  the  solid  line  Is  the  prediction  of  the  model,  and  the 
circles  are  experimental  points  from  some  BRL  «rark  [26] . 
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